Shear thickening and jamming in densely packed suspensions of different particle 

shapes 
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We investigated the effects of particle shape on shear thickening in densely packed suspensions. 
Rods of different aspect ratios and non-convex hooked rods were fabricated. Viscosity curves and 
normal stresses were measured using a rheometer for a wide range of packing fractions for each 
shape. Suspensions of each shape exhibit qualitatively similar Discontinuous Shear Thickening. 
The logarithmic slope of the stress/shear-rate relation increases dramatically with packing fraction 
and diverges at a critical packing fraction (j)c which depends on particle shape. The packing fraction 
dependence of the viscosity curves for different convex shapes can be collapsed when the packing 
fraction is normalized by (j)c. Intriguingly, viscosity curves for non-convex particles do not collapse 
on the same set as convex particles, showing strong shear thickening over a wider range of packing 
fraction. The value of (j)c is found to coincide with the onset of a yield stress at the jamming 
transition, suggesting the jamming transition also controls shear thickening. The yield stress is 
found to correspond with trapped air in the suspensions, and the scale of the stress can be attributed 
to interfacial tension forces which dramatically increase above (j)c due to the geometric constraints 
of jamming. Using this connection we show that the jamming transition can be identified by simply 
looking at the surface of suspensions. The relationship between shear and normal stresses is found to 
be linear in both the shear thickening and jammed regimes, indicating that the shear stresses come 
from friction. In the limit of zero shear rate, normal stresses pull the rheometer plates together due 
to the surface tension of the liquid below (/)c, but push the rheometer plates apart due to jamming 
above (j)c. 

PACS numbers: 83.80.Hj, 83.85.Cg, 83.60.Fg 



I. INTRODUCTION 

Shear thickening is a category of non-Newtonian fluid 
behavior in which the viscosity increases as a function 
of shear rate. A particularly dramatic form, known as 
Discontinuous Shear Thickening, occurs in many densely 
packed suspensions and colloids, the most well-known ex- 
ample being cornstarch in water. These suspensions feel 
like a thin liquid when stirred weakly, but feel very thick 
when stirred harder, and feel thin again when the stress 
is removed. Discontinuous Shear Thickening is quantita- 
tively characterized by a steep jump in the shear stress 
r with increasing shear rate 7, which becomes steeper 
with increasing packing fraction up to the point where it 
appears to be discontinuous The stress jump can 

be attributed to frictional dissipation from chains of solid 
particle contacts that span the system due to boundary 
confinement in response to dilation under shear [5] . Here 
we investigate how Discontinuous Shear Thickening de- 
pends on the shape of particles in suspension. Designer 
particle suspensions are of potential practical interest due 
to the damping and shock absorbing abilities of Discon- 
tinuous Shear Thickening suspensions [6]. 

Different particle shapes also provide an opportunity 
to investigate the relation between jamming and Discon- 
tinuous Shear Thickening, since different particle shapes 



can pack at very different volume fractions ^[7,3]- While 
it has long been suggested that shear thickening is related 
to jamming [9], there have been few quantitative connec- 
tions. In a previous letter, we showed that the packing 
fraction dependence of the slope of r(7) can be char- 
acterized as a power law diverging at a critical packing 
fraction (j)c for nearly-spherical particles [10]. Above (pc 
these systems are jammed, meaning they will not flow for 
applied stresses below a yield stress because geometric 
constraints from particle contacts suppress the relative 
shearing of particles [11^. In this paper we extend our 
work in Ref. [10] to further characterize the relationship 
between the mechanics of shear thickening and jammed 
systems using different particle shapes. 

A few previous works have investigated effects of par- 
ticle shape on shear thickening. Clarke [l2| compared 
shear thickening for different shapes, but this was done 
at a single packing fraction and so did not relate shear 
thickening to the packing behavior of the particles. Bea- 
zley [13] showed that higher aspect ratio particles shear 
thickened at lower packing fractions, and suggested that 
this was related to the lower random packing densities 
of the higher aspect ratio particles. More recently, Egres 
& Wagner [8] performed detailed measurements of shear 
thickening for colloidal rods of different aspect ratio. 
They compared the packing fraction range of Discontinu- 
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ous Shear Thickening with dense packings obtained from 
centrifugation, and observed a similar qualitative trend of 
decreasing packing fraction with aspect ratio. Here we go 
beyond those works in two important ways. First, using 
the unique PRINT® process [14], we were able to fabri- 
cate out of identical materials a wider variety of particle 
shapes including rods with different aspect ratios, as well 
as some non-convex rods with hooks. Second, we show a 
quantitative match between critical points for shear thick- 
ening and jamming by quantifying the slopes of the shear 
thickening regime, the yield stress, and normal stresses 
as a function of packing fraction for different shapes. 

The rest of this paper is organized as follows. Rhe- 
ological measurement methods and particle details are 
described in Sec. [Ill Viscosity curves for particles of dif- 
ferent shapes over a range of packing fractions are shown 
in Sec. Illli We then compare the strength of shear thick- 
ening for different shapes based on the logarithmic slopes 
of r(7) as a function of packing fraction cf) for differ- 
ent particle shapes in Sec. IIII A\ where we show that 
the packing fraction dependence for many convex shapes 
can be collapsed relative to the shape-dependent criti- 
cal packing fraction (pc- In Sec. IIV Al we show that yield 
stresses which are mainly measured above 0c are corre- 
lated strongly with trapped air in the suspensions sug- 
gesting the stresses come from interfacial tension. In 
Sec. IIVBI we use images of the surface of suspensions 
at different packing fractions to further demonstrate the 
connection between surface tension and the yield stress. 
In Sec. |Vl we show measurements of normal stress in re- 
lation to shear stress suggesting a frictional relation, and 
how the static normal stress changes abruptly at 0c- In 
Sec. lVII we use the measurements of slopes of the viscosity 
curves for different packing fractions and shapes to rein- 
terpret previous results that were attributed to particle 
shape or polydispersity. 



II. MATERIALS AND METHODS 



/im) by (24.9 =b 0.6 fim) by (30.2 ± 0.4 fim) referred to 
as aspect ratio F = 6 rods, and the third has dimensions 
(21.3 ± 2.5 /im) by (25.3 ± 0.8 /im) by (25.3 ± 0.8 /im) 
referred to as aspect ratio F = 1 rods. A fourth shape 
was fabricated, shown in Fig. [T]i, which is a variation on 
the F = 9 rods with cubic "hooks" on opposite ends, so 
the particle has the shape of an 'S'. These are referred to 
as hooked rods. These rods have dimensions (198.1 ± 1.1 
/im) by (20.6 ±0.4 /im) by (24.5 ±0.8 /im), and the cubic 
hooks are 23.0 ±1.2 /im long, share the width (20.6 ± 0.4 
/im) with the rod, and are 24.7 ± 1.2 /im wide along the 
length of the rod. 

The particles were suspended in poly(ethylene glycol) 
dimethyl ether (PEG) (Mn = 250g/mol) with a viscos- 
ity of r]i = 60 mPa s. The particles are nearly density 
matched with the PEG solvent; based on settling times 
we estimate the density of the particles to be 1% greater 
than that of the fluid [16]. 

For a comparison to standard spherical particles, and 
for some direct visualizations, we used opaque polyethy- 
lene spheres obtained from Cospheric. These particles 
have a nominal diameter range of 125-150 /im and den- 
sity of 1.01 g/mL. They were dispersed in silicone oil AR 
20 with a nominal density of 1.01 g/mL and viscosity 
of 20 mPa s. This suspension is then density matched 
±0.01 g/mL. 

B. Methods 

Viscosity measurements were made with an Anton 
Paar Physica MCR 301 rheometer. The sample of thick- 
ness d is held in a thin layer between two horizontal cir- 
cular smooth plates of radius R = 12.5 mm. The sides 
of the sample are held at the edge of the plates by sur- 
face tension. The rheometer applies a constant torque 
T which shears a sample while recording the top plate 
rotation rate uj. The viscosity is measured as 7^ = r/7 in 
a steady state for shear stress 



A. Materials 

For this study we fabricated particles on the scale of 20- 
200 /im to ensure geometric effects were prominent over 
Brownian motion and interparticle interactions. We also 
designed the chemistry specifically to minimize particle- 
fluid surface tension so the suspension would exhibit 
shear thickening [15]. The particles were fabricated us- 
ing the PRINT® process [14], in which molds are made 
to design the shape of particles. The particles were com- 
posed of 57% (w/w) triacrylate, 40% (w/w) methyl ether 
acrylate, 1% (w/w) 1-hydroxycyclohexyl phenyl ketone, 
and 2% (w/w) fluorescein o-acrylate. Details of the pro- 
cess used to make the particles for this study can be found 
in Ref. [l6[. Three types of rectangular rods were used, 
shown in Fig.[TJ The first has dimensions (266.5±0.6 /im) 
by (29.9±0.5 /im) by (31.4±0.4 /im) referred to as aspect 
ratio F = 9 rods, the second has dimensions (139.5 ± 0.7 
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The normal stress tn is measured as the normal force 
pushing against the top plate divided by the cross- 
sectional area of the plate. These definitions are meant 
to characterize the average mechanical response due to 
an external force in a size-independent way. In systems 
that exhibit Discontinuous Shear Thickening, this defini- 
tion can differ dramatically from the local hydrodynamic 
constitutive relation between shear stress and shear rate 
because the shear stress is strongly coupled to normal 
stresses which depend on boundary conditions ^] . 




FIG. 1: Scanning electron microscope images of dry PEG particles, 
rods, (d) r = 9 hooked rods 



(a) Aspect ratio T = 9 rods, (b) T = 6 rods, (c) T = 1 



Measurements were made at a bottom plate tempera- 
ture controlled at 20° C with the room humidity rang- 
ing from 22% to 38%, although during individual exper- 
iments the humidity was constant. This affects solvent 
evaporation/adsorption which can have a significant ef- 
fect on the rheology due to the sensitive packing fraction 
dependence of shear thickening suspensions [lo| . To min- 
imize evaporation or adsorption, we used a solvent trap 
which enclosed the sample and a small amount of air 
around it by an extra layer of liquid. 

Packing fractions (j) were defined as the volume of solid 
particles over the total volume of solids plus liquid. Each 
ingredient was measured by mass which was divided by 
the material density to obtain the volume. The particles 
and liquids were stirred together for 30 s using a vortex 
mixer. The packing fraction in terms of the inverse of 
the available free volume per particle decreases slightly 
during measurements as the samples are seen to dilate 
during shear thickening [5|. There may also be a differ- 
ence in packing fraction definitions due to air bubbles 
becoming trapped in the interior of the suspension (see 
Sec. IIV A|) . Because of the small size of our samples and 
changes over time due to humidity, our absolute measure- 
ments of packing fraction could vary by as much as 0.04 
for each shape upon trying to reproduce mixtures. This is 
the relevant error for comparing to external experiments. 
Within an individual measurement series where the pack- 
ing fraction was lowered by adding a small amount of liq- 
uid PEG for each successive experiment, the uncertainty 
in packing fraction when comparing to neighboring data 
points within that series was estimated to be 0.005. 

Samples were pre-sheared immediately before measure- 



ments for at least 100 seconds at shear rates above the 
shear thickening regime where the steady state flow is 
fully mobilized [5]. After this pre-shear, measurements 
on suspensions were found to be repeatable with a typical 
variation of 10-20% from run to run. This is the variation 
whether we remeasure a sample that is in place on the 
rheometer or replace it with a new one with the same 
procedure. Following the preshear, measurements were 
performed with slowly decreasing and then increasing 
stress ramps and additional runs were made with differ- 
ent control ramp rates to check for hysteresis, thixotropy, 
and transients. Examples of hysteresis loops for different 
ramp rates are shown in Ref. In this paper, viscos- 
ity curves are shown for a ramp rate of 500 s per decade 
of stress which is slow enough to be in the steady state 
limit. We show only one set of curves for brevity since 
in the steady state limit they were all identical within 
typical variations. 



Because only small samples of particles could be fab- 
ricated, the sample thickness d was limited to no more 
than a few particles across. We ensured that in each 
experiment, d was more than twice the longest length 
the particle, which is large enough to avoid significant 
finite-size effects on the viscosity curves [16]. Slip was 
shown not to have a significant effect on viscosity curves 
for similar measurements with some of the same parti- 
cles [16]. The suspensions were observed not to spill out 
of the rheometer until shear rates exceeding those of the 
reported measurements; however in most cases this set 
the upper limit of the range of measurements. 
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FIG. 2: (color online) Viscosity vs. shear stress for rods of 
aspect ratio F = 9. Packing fractions are decreasing from 
top to bottom. The solid line has a slope of 1 corresponding 
to a constant shear rate. The dashed line in the key indicates 
the critical packing fraction (\)c above which the suspension is 
jammed with a large yield stress, and below which the sus- 
pension exhibits shear thickening. 
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FIG. 3: (color online) Viscosity vs. shear stress for rods of 
aspect ratio F = 6. Packing fractions are decreasing from 
top to bottom. The solid line has a slope of 1 corresponding 
to a constant shear rate. The dashed line in the key indicates 
the critical packing fraction (\)c- 



III. VISCOSITY CURVES 

In this section we show viscosity curves as viscosity 
vs. shear stress for suspensions of particles with dif- 
ferent shapes over a range of packing fractions. Vis- 
cosity curves for rods of aspect ratios F = 9, F = 6, 
F = 1, and F = 9 hooked rods are shown in Figs. [21 [3l 
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FIG. 4: (color online) Viscosity vs. shear stress for rods of 
aspect ratio F = 1. Packing fractions are decreasing from 
top to bottom. The solid line has a slope of 1 corresponding 
to a constant shear rate. The dashed line in the key indicates 
the critical packing fraction (\)c. 
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FIG. 5: (color online) Viscosity vs. shear stress for F = 9 
hooked rods. Packing fractions are decreasing from top to 
bottom. The solid line has a slope of 1 corresponding to a 
constant shear rate. The dashed line in the key indicates the 
critical packing fraction (\)c. 



[31 and [5] respectively, for different packing fractions ^. 
In each case, the viscosity curves exhibit a shear thick- 
ening regime indicated by a positive slope in some in- 
termediate stress range which does not vary too much 
with packing fraction. The slope of this shear thickening 
regime becomes steeper with increasing packing fraction, 
ultimately reaching a maximum slope of about 1 on a log- 
log plot of viscosity vs. stress. Because of the definition 
7^ = r/7 this limit corresponds to a constant shear rate 




5 



or equivalent ly an infinite slope of shear stress vs. shear 
rate. At higher packing fractions, instead of shear thick- 
ening, a large yield stress is observed, which is indicated 
by the viscosity diverging as the shear stress decreases. 
This behavior is typical of Discontinuous Shear Thicken- 
ing in nearly-spherical particles [1, HI, [lo| or rods @ . 

One notable feature of the viscosity curves of Figs. [21 
[31 [H and [5] that is dependent on particle shape is the 
variation in the packing fraction (j)c where jamming oc- 
curs. Jamming is usually characterized by the onset of 
a yield stress in ideal systems [18]. Practically, we mea- 
sure the onset of jamming as the packing fraction where 
the largest increase in the yield stress occurs [10], indi- 
cated in by the dashed lines in the keys of Figs. [2 [3l 
[H and [5l The very small yield stresses at lower packing 
fractions are too weak to attribute to geometric jamming 
(see Sec. [IV] for a detailed investigation). While suspen- 
sions of frictional spheres jam at about (pc = 0.58 [1(|, 
we find = 0-55 for F = 1 rods, <pc = 0-37 for F = 6 
rods, (j)c = 0.35 for F = 9 rods, and (j)c = 0.34 for the 
hooked rods. The trend of decreasing packing fraction 
with increasing rod aspect ratio is consistent with other 
measurements 0, 111. 

A. Packing fraction dependence of slopes 

To quantify the strength of shear thickening, we con- 
sider the steepness of the viscosity curves following our 
method in Ref. ^]. For each viscosity curve, we fit 
T] (X T^~^ (equivalent to r ex 7^^^) to the shear thicken- 
ing regime. The value of e characterizes the strength of 
shear thickening such that the bound e = 1 corresponds 
to a Newtonian fluid and the bound e = corresponds 
to the steepest possible shear thickening curve with an 
infinite slope of stress vs. shear rate. In Ref. [lo| we 
measured that e goes to zero in the limit as (j) approaches 
(pc from below for both suspensions of glass spheres and 
cornstarch particles. The critical packing fraction 0c was 
also found to occur at the same packing fraction as the 
divergence of the magnitude of the viscosity and the on- 
set of the yield stress. In cases where there is a large yield 
stress that hides much of the shear thickening regime so 
the stress-shear rate relation did not exhibit a constant 
slope over a wide range, we had to modify this procedure 
to measure the underlying differential increase in stress 
with shear rate. In such cases we fit 

^(7) =ry^ ai7^/2 + a27'/' (3) 

to the data in the steep part of the shear thickening 
regime and at lower stresses [15]. The first two terms 
of Eq. [3] characterize the yield stress and shear thinning 
behavior as in a Herschel-Bulkley model. When the yield 
stress and correspondingly the first two terms of Eq. [3] 
are small, this procedure is equivalent to the original 
procedure of fitting just the shear thickening regime to 




0.4 0.6 0.8 1.0 

FIG. 6: (color online) The strength of shear thickening is 
characterized by a fit of r oc 7^^^ to the shear thickening 
regime, and e is plotted vs. normalized packing fraction 0/0c 
for different particle shapes. Black solid squares: rods with 
aspect ratio F = 9, (/)c = 0.35. Red diamonds: rods with 
aspect ratio F = 6, (/)c = 0.37. Green open squares: rods 
with F = 1, (j)c — 0.55. Blue crossed squares: hooked 
rods, (j)c — 0.34. Black open circles: glass spheres in water, 
(j)c — 0.58 [13] • Purple down-pointing triangles: cornstarch in 
glycerol, 0c = 0.57 [^. Orange up-pointing triangles: corn- 
starch in water, (/)c = 0.48 ^]. The collapse of the data in 
the gray band suggests that the normalized packing fraction 
(j)/(t)c determines the strength of shear thickening for convex 
particle shapes. The dashed line at e = 1 corresponds to vis- 
cous flow with T oc 7, the limit at small (j) if inert ial effects 
are negligible. The dashed-dotted is e = 2/3, corresponding 
to an inertial flow with r oc 7^^^. The solid line is a best fit of 
e oc {(j)c — (j))^ to the data for convex shapes and (j)/(j)c > 0.8. 



Here we apply this same analysis to suspensions of 
particles of different shapes. In Fig. [6] we plot e vs. a 
normalized packing fraction (j)/(j)c' The value of (j)c is 
measured independently for each shape as the packing 
fraction where there is a steep jump in the yield stress. 
Data is shown for each of the shapes studied here, as 
well as glass spheres from Ref. [10] and cornstarch from 
Ref. [5*]. With the exception of the hooked rods, the data 
for every convex shape collapse near highlighted by 
the shaded band in Fig.[6l This suggests that the normal- 
ized packing fraction 0/0c generally characterizes shear 
thickening for each of those convex shapes, despite the 
large differences in 0c- Since the critical packing frac- 
tion 0c is measured based on the yield stress, which is 
independent from the slope of the viscosity curve, the 
collapse confirms a connection between shear thickening 
and the jamming transition. Specifically, the fact that 
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the slope diverges (e ^ 0) as goes to (/)c suggests the 
proximity to the jamming transition controls the strength 
of shear thickening like a critical point in a second order 
phase transition. The scatter from repeated runs is simi- 
lar to the scatter around zero near (/)c, indicating the typ- 
ical uncertainty in fitting the viscosity curves to obtain 
e. While there has been some success in characterizing 
the order parameter and scaling exponents for the jam- 
ming transition j^, how critical scaling theory applies 
to Discontinuous Shear Thickening is still an open prob- 
lem. The situation is complicated by the fact that these 
suspensions are neither in equilibrium nor homogeneous, 
nor does the divergent scaling apply to local constitutive 
laws [El. 

The collapse of the slopes in Fig. [6] suggests that some 
dense suspensions that show strong shear thickening but 
do not exhibit a discontinuous stress/shear-rate relation 
(e 0) as in Fig. 3] may still fit into the same model as 
Discontinous Shear Thickening suspensions. It has been 
shown that the stress contribution from the yield stress 
can be added linearly to the stress contribution respon- 
sible for Discontinuous Shear Thickening [15]. Since the 
total viscosity from the sum of these two terms is defined 
as an absolute ratio r] = r/j, the yield stress can make a 
significant contribution at higher shear rates and the to- 
tal viscosity can be less strongly shear thickening than the 
differential increase in stress with shear rate from shear 
thickening mechanisms. Based on this interpretation, the 
large yield stresses near (j)c in Fig.|4]hide much of the oth- 
erwise shear thickening regime, and the apparent slope 
of the viscosity curve is lessened by a broad crossover 
from the dominance of the yield stress to the dominance 
of the shear thickening term. This could mean that the 
stresses from shear thickening mechanisms are the same 
as in other dense suspensions, but if the stresses from 
yield stress mechanisms are large they can be dominant 
in the overall rheology. 

To quantify the possible critical scaling behavior, we fit 
the function e oc (0* — to all of the data for the convex 
shapes in Fig. [6] for ^/^c > 0.8. Statistical uncertainties 
of 0.015 on ^ and 0.09 on e are based on repeated mea- 
surements. Varying the fit value of the exponent ^, 
and the proportionality, we obtain ^*/^c = 0.979 ±0.004 
and ^ = 0.52 ± 0.06. Limiting the range of the fit to 
higher ^, or varying the relative uncertainties on pack- 
ing fraction and e, does not change the fit values signif- 
icantly. Because we can only obtain measurements of e 
when ^ < only negative errors in (/) are measured 
near which biases the distribution of data near the 
critical point so the best fit value could be smaller than 
the actual value by about the uncertainty on (j). Given 
this bias, the packing fraction where the viscosity curve 
becomes discontinuous and e goes to zero (represented 
by ^*) is consistent with the packing fraction where the 
large yield stress onsets. If we then redo the fit while 
fixing 0* = 0c, then we obtain ^ = 0.74 ± 0.09. 

Further away from 0c , the critical point no longer dom- 
inates the rheology and different forces can be dominant. 
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FIG. 7: (color online) (a) Scatter plot of yield stress normal- 
ized by the surface tension stress scale 7/a vs. normalized 
packing fraction (/)/(/)c for PEG particles with different shapes 
identified in the key, along with data for polyethylene spheres 
(pink circles), (b) Normalized difference between the volume 
of the suspension in the rheometer gap Vg and the volume 
of solids and liquids mixed in suspension Vs, indicating the 
fractional volume of air trapped in suspension. Dashed line: 
minimum fraction of trapped air due to geometric constraints 
for hard particles. The absolute uncertainty in the volume 
ratio on the vertical scale is 0.02 except for the spheres. 



For higher viscosity fluids, viscous forces can dominate, 
where a Newtonian behavior is expected with e = 1 
(shown by the dashed hue in Fig. [6j), which is satisfied 
for example by cornstarch in a glycerol-water mixture 
at small (j) {r]i = 80 mPa s, data from Ref. Q). For 
lower viscosity fluids, inert ial forces can dominate in the 
stress regime of our measurements. For cornstarch in 
water {rji = 1 mPa s), this results in an intermediate 
Reynolds number regime in the stress range of interest 
where r ex 7^/^, corresponding to e = 2/3 (shown by the 
dashed-dotted line in Fig. [6l data from Ref. [5]). 



IV. YIELD STRESS 

Here we address the origin of the yield stress seen in 
the viscosity curves in Figs. [21 [3l HI and O We show 
a scatter plot of yield stress Ty vs. normalized packing 
fraction 0/(/)c in Fig. [7^. The large yield stress above 0c 
suggests a jammed state, which is attributed to geometric 
constraints on particle motion when the packing fraction 
becomes so high that particles cannot move around each 
other In systems of confined volume, the scale of 

the yield stress is traditionally attributed to the stiffness 
of particles |l8|. However, for our suspensions, the mea- 
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sured value of the yield stress is only enoug h to co mpress 
hard elastic spheres by a strain of order ^JryjE^ < 10~^ 
using a Hertzian contact model where is the Young's 
modulus of the particles [18]. Under this model, our sus- 
pensions could not be compressed to packing fractions 
above the jamming transition by more than 10~^. The 
smaller yield stresses for (/) < in Fig. [7^ also lacks an 
explanation. 



A. Trapped air 

The key difference from those models is that our sus- 
pensions are not confined to a fixed volume; instead they 
are confined by surface tension at the liquid-air interface 
which is much softer than the particles. If the particles 
are packed to a density just above the jamming transi- 
tion, then the liquid-air interface could deform by a dis- 
tance up to about half a characteristic particle diameter 
a (calculated as the cube root of the volume) to accom- 
modate the particle packing. This finite-size effect would 
only allow the measured packing fraction to reach up to 
only ajR ~ 0.004 above jamming in our case. These 
limitations imply that even if we try to mix higher ra- 
tios of solids to liquid plus solids, there is not enough 
confining stress to pack the particles much more densely 
than at the jamming transition, and so the particle pack- 
ing must take up more volume by trapping air bubbles 
within the packing. We propose that the interfacial ten- 
sion from these air bubbles and the deformed liquid-air 
interface at the boundary can be a source of yield stress. 
A curved liquid-air interface produces a force from inter- 
facial tension on the attached particles. If solid frictional 
contacts between particles span the system, these forces 
can transmit along chains and be redirected in a ran- 
dom packing such that the resistance to shear stress is 
proportional to the normal stress due to confinement as 
in a frictional scaling. Such force chains are expected 
at packing fractions above the jamming transition [9], at 
lower packing fractions with even tiny interaction forces 
between particles [l9| , and during shear thickening when 
packings dilate [sj. The scale of the stress from surface 
tension scales as 7/r for a surface with curvature of ra- 
dius r. Presuming the scale of the curvature is set by the 
particle size, this gives a stress scale 7/a which we can 
compare yield stress measurements to. 

We measure the volume of trapped air based on the 
difference between the volume that the suspension takes 
up in the rheometer gap Vg with an uncertainty of 2%, 
and the volume of solids plus liquids mixed in suspen- 
sion Vs. We show this difference in Fig. [TId as a function 
of packing fraction for different particle shapes, with the 
difference in volumes normalized as V^/Fs — 1 indicating 
the fractional volume of air over liquids plus solids in the 
suspension. It can be seen that the amount of trapped air 
is consistent with zero at low packing fractions, but can 
become somewhat positive as approaches 1 from 

below and typically jumps to large positive values for 



^/0c > 1- This confirms that there is air trapped in 
the suspensions above jamming. The amount is consid- 
erably more than the minimum required by geometric 
constraints on hard particles which will not pack more 
densely than at the jamming transition (dashed line in 
Fig. [7]). Additionally, it is seen that some air tends to 
become trapped in the suspensions even below jamming 
but at high packing fractions where the effective viscosity 
happens to be large. The tendency for air to be trapped 
in suspensions that are jammed or nearly so is likely due 
to the difficulty of bubbles escaping between particles 
that are closely packed in suspensions with a high vis- 
cosity or yield stress. We also plot in Fig. [71 the yield 
stress and fraction of trapped air for suspensions of 135 
/im polyethylene spheres in mineral oil whose volumes 
were measured much more precisely. These samples were 
sonicated after mixing to increase the mobilization of par- 
ticles and help trapped air bubbles escape. These data 
are seen to match more closely to the minimum amount 
of trapped air plotted as the dashed line, which suggests 
the amount of trapped air can vary with the mixing pro- 
cedure. 

By comparing panels a and b of Fig. [71 it can be 
seen that there is a strong correlation among the dif- 
ferent shapes between the yield stress and the amount 
of trapped air for (j) < (j)c. In particular, the particle 
with both the largest yield stress and largest amount of 
trapped air, the F = 1 rods, is the smallest particle used, 
which might enhance the tendency to trap air bubbles 
because stresses from surface tension are relatively larger 
for smaller particles. 

To connect the trapped air to the yield stress, we show 
a scatter plot of yield stress Ty vs. trapped air fraction 
Vg/Vs — 1 in Fig. [8] The data correspond to different 
packing fractions and particle shapes. The yield stress 
values are normalized by the predicted interfacial stress 
scale 7/a. A strong positive correlation is seen between 
the yield stress and volume of trapped air, and the data 
for different shapes even collapse in a narrow band. The 
data clustered at the lower left generally correspond to 
(j) < (j)c where the yield stress is small or below the reso- 
lution limit and there is little or no measurable trapped 
air. The data above 0c where the yield stress is large 
generally have large volumes of trapped air so show up 
on the upper right of the figure. The scale of the yield 
stress in the limit of larger volumes of trapped air is rem- 
iniscent of other systems jammed by interfacial tension. 
In jammed foams, the yield stress is observed to be on 
the scale of Ty ~ 0.057/a [20]. In 3-phase suspensions 
with 2 non-miscible fluid phases, a yield stress was ob- 
served to be on the scale of Ty ^ 0.057/a for moderate 
fractions of the 3rd phase [21]. While the 3rd phase in 
that case was a liquid instead of gas, the forces from in- 
terfacial tension are much stronger than gravity on the 
scale of the particle size. The maximum stress seen in the 
shear thickening regime is also found to be on the order 
of 0.057/a [5]. While in the latter case the suspensions 
are not statically jammed, the measured forces are simi- 
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FIG. 8: (color online) Scatter plot of yield stress vs. volume 
of air trapped in rheometer as a fraction of solid plus liquid 
sample volume for a variety of PEG particle shapes listed in 
the key. Data for polyethylene spheres (pink circles). For a 
given series of one particle shape, packing fraction generally 
increases to the right, and the cluster of points in the lower 
left corresponds to packing fractions below jamming. Dashed 
line: 0.057/a corresponding to the limiting stress value seen 
in foams jammed by interfacial tension [20j], 3-phase suspen- 
sions oi" the maximum stress seen in the shear thickening 
regime Q]. 



lar because they are transmitted through the suspension 
along solid contacts between particles, which form when 
the suspension dilates under shear against the interfacial 
tension at the boundary. The similarity with these other 
results strongly suggests that the value of the yield stress 
in these jammed suspensions comes from interfacial ten- 
sion. 

While the specific values of the yield stress and vol- 
ume of trapped air shown in Fig. [7] likely depend on 
the sample preparation procedure, the collapse of yield 
stress data for the PEG particles in Fig. [8] suggests a 
general relation between the yield stress and fraction of 
trapped air, despite the different particle shapes. If the 
yield stress is given by the average capillary stress 7/r 
percolated across the suspension along force chains, then 
the data collapse suggests the average ratio of particle 
size to bubble size is similar from shape to shape. While 
it is not obvious whether the characteristic bubble size is 
determined by hydrostatics or shear, we note that one of 
the contributions of a preshear in bringing a sample to 
a reproducible state could be in bringing the bubble size 
distribution to a steady state. 



B. Surface roughness and capillary forces 

If a suspension is jammed, the boundaries must be able 
to balance the forces that are transmitted along force 
chains through the bulk. When the boundary is a liquid 
air interface, it must then become curved so that surface 
tension can provide the stress to push back on particles 
which penetrate the liquid-air interface . If the parti- 
cles are between about 1 and 100 /im, then the grains are 
large enough to scatter light diffusively and small enough 
that they cannot be seen individually, so the surface ap- 
pears rough by eye. On the other hand, if a suspension 
is unjammed there is free space for the particles to rear- 
range, so any particles on the surface will be pushed by 
surface tension to the interior (assuming the liquid wets 
the particles, which is also a requirement to observe shear 
thickening [15|), resulting in a flat and shiny surface. In 
this section we demonstrate this contrast in surface ap- 
pearance between jammed an unjammed states and that 
it can be used to make quantitative measurements of the 
jamming transition ^c- 

The suspensions used consisted of 135 /im polyethylene 
spheres density matched in silicone oil, which are opaque 
and do not settle. These suspensions were placed in a 
H = 1.3 mm deep layer and viewed from the top with di- 
rect lighting as shown in Fig. [9l The series of images was 
taken by starting at a high packing fraction well above 
the jamming transition (0 > 0c ), then adding oil to re- 
duce the packing fraction for each successive image. At 
each packing fraction the sample was vibrated at 40kHz 
with a sonicator for 1 min, tilted during and after sonica- 
tion to observe whether or not shear occurred, then any 
flow was allowed to come to rest before the image was 
taken. 

For the highest packing fractions (j) > 0.64, the sus- 
pension surfaces appear dry and rough and did not flow 
under tilt during or after sonication. For 0.64 > > 0.57 
they appear wetter but still rough. In this range the sam- 
ples did flow while being sonicated, but afterwards did 
not flow at inflnitesimal tilt angles (measured at ~ 1°). 
However, for several of the packing fractions in this range 
there was a critical tilt angle O above which flow could 
be found. For (j) < 0.57, the samples appeared very shiny 
and smooth, and they flowed even at inflnitesimal tilt 
angles. 

The critical packing fractions for the transitions in 
Fig. [9] match up with those usually found for the jam- 
ming transition in similar systems [ll|. The transition 
at (/) = 0.64 corresponds to random close packing for 
frictionless spheres, above which packings of spheres are 
jammed and below which they can flow [22]. When com- 
paring to other published values for critical packing frac- 
tions, it is reasonable to expect an absolute uncertainty 
in the range of 0.01 in the packing fraction due to factors 
such as a flnite size effect and polydispersity in the par- 
ticle sizes, and sample preparation. Experimentally it is 
usually found that packings of spheres remain mechani- 
cally stable, i.e. have a yield stress down to packing frac- 
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FIG. 9: (color online) Top views of a static layer of 135 fxm polyethylene spheres density matched in silicone oil. The layer was 
1.3 mm deep and each image is 2 cm on a side. Images were taken for different packing fractions given in the corner of each 
panel. Lighting was direct to emphasize changes in reflectivity. Packing fractions > 0.64 did not flow under any amount of 
tilt. Packing fractions 0.64 > (j) > 0.58 flow while being sonicated, but after stopping sonication, these suspensions no longer 
flowed at small tilt angles. Packing fractions < 0.57 flow under even small tilt 1 degree). The reflectance of the surface 
increased so much when was decreased to 0.57 that the camera became oversaturated by the direct reflection of the light. 



tions as low as 0.56, called random loose packing, with 
the value of the packing fraction depending on friction 
and the density difference between the particles and fluid 
|25|, l26j. Packing fractions closer to (j) = 0.64 can only 
be reached for frictional, settling particles if the packings 
are vibrated which mobilizes the particle contacts, effec- 
tively eliminating the effect of friction [23^ . This explains 
why in the range 0.64 > <j) > 0.57 our suspensions flow 
under sonication but are otherwise jammed. This easily 
visible indicator of jamming based on surface roughness 
contrasts with other jammed systems, where there is no 
such visual test. The insterstitial liquid allows this sen- 
sitive measurement by eye because changes in the fluid 
level at the surface due to changes in packing density only 
have to be on the scale of a particle size to dramatically 
change the surface appearance. 

We note the visible transition at = 0.64 in Fig. [9] 
is not always observed, depending on sample prepara- 
tion. If we do not sonicate, the samples can appear to be 
dry until they are diluted all the way down to ^ = 0.57. 
This could be because without the mobilization of par- 
ticle contacts allowed by sonication, the suspension may 
not be able to pack more densely than (j) = 0.57, perhaps 
trapping air bubbles, or to make up the extra space the 
liquid may retreat to the interior of the sample. Despite 
this hysteresis effect in the visible transition at <p = 0.64 
during preparation, rheology measurements with pres- 
hear consistently show the onset of a yield stress at the 
same packing fraction. 

We can connect the value of the critical tilt angle O for 
flow to the yield stress of the suspension. For a tilt angle 
0, the stress applied by gravity parallel to the surface is 
pgH sin O, which is shown in Fig.[TO]for packing fractions 
in which one was measured. We compare this to the yield 
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FIG. 10: Open symbols: gravitational stress pgH sin S re- 
quired to initiate shear under tilt for the suspensions shown 
m Fig.H The minimum tilt angle G required for shear is mea- 
sured relative to a horizontal surface. Dashed line: threshold 
packing fraction above which the sample did not flow at any 
tilt angle (without sonication). Solid symbols: yield stress Ty 
obtained from viscosity curves in the rheometer. The jam- 
ming transition at (j) = 0.57 above which there is a yield 
stress coincides with the visible change in the surface shown 
in Fig. El 



stress Ty obtained from rheometer measurements in the 
zero shear rate limit. For (j) < 0.58, we measured no 
yield stress, within our experimental resolution of 10~^ 
Pa. There is an absolute uncertainty of about 0.005 in 
the packing fraction measurements for the yield stress 
due to the process of loading the sample. Within this 
uncertainty, the shear stress provided by gravity at the 
onset of shear under tilt matched with the measured yield 
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FIG. 11: (color online) Normal stress tn vs. shear stress r for 
r = 9 rods. Packing fractions (j) are shown in the key. Solid 
line: linear relationship corresponding to a frictional response. 
The dashed line in the key indicates the critical packing frac- 
tion (/)c above which the suspension is jammed with a large 
yield stress, and below which the suspension exhibits shear 
thickening (see Fig. [2|). 



stress in the rheometer. This confirms that the visible 
changes in the surface can be quantitatively matched to 
changes in the yield stress. 

The above analysis implies the macroscopic roughness 
of the surface of a suspension can be used as an indica- 
tor of the yield stress. At each packing fraction shown 
in Fig. [9]for ^ > 0.57, the sample surface has roughness 
on a macroscopic scale, i.e. much larger than individual 
particles. If an asperity of height H forms on the upper 
surface of a fluid it can remain stable if the yield stress 
exceeds a value on the scale of pgH. For these samples 
we observed asperities on the order of 1 mm at the high- 
est packing fractions, consistent with the measured yield 
stresses on the order of 10 Pa and the value H used in 
the gravitational stress scale. 

The indication of a yield stress by the rough appear- 
ance of the surface confirms that surface tension forces 
are confining these suspensions at the boundary, whether 
or not there is trapped air. Since the scale of forces from 
surface tension are similar at the boundary or due to 
trapped air, in general it could be possible that the de- 
formation at the liquid-air interface at the boundary is 
responsible for the yield stress, which is then transmit- 
ted along solid particle contacts through the suspension. 
For our measurements, we know there is also trapped air 
based on volume measurements so the deformation of the 
liquid-air interface at the boundary is also balancing the 
stresses from surface tension in the interior. 
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FIG. 12: (color online) Inverse friction coefficient from 
a fit of Tn — tn,q + iJi~^T. Data are shown for several parti- 
cle shapes listed in the key vs. normalized packing fraction. 
Significant values of indicating a frictional relation, are 
observed for cj) ^ 0.8(/)c which is the same range where Discon- 
tinuous Shear Thickening is seen in Fig. [G] 

V. NORMAL STRESSES 

In this section we show the relation between the normal 
stresses and shear stresses. We measured normal stresses 
along with the viscosity curves of Figs. El [3l HI and [5l 
Some examples of normal stresses are plotted vs. shear 
stresses in Fig. [11] for the F = 9 rods at different pack- 
ing fractions. For packing fractions below (/)c, the normal 
stress at zero shear stress is negative (pulling on the top 
plate), while above the normal stress is everywhere 
positive (pushing on the top plate). For each packing 
fraction, we see an increase in normal stress with increas- 
ing shear stress. The roughly linear slopes and the fact 
that the shear and normal stress are on the same order 
of magnitude suggests a frictional scaHng. 

To quantify such a frictional scaling, we use a fit- 
ting function tn = tn^o + /i~^T with fit parameters 
/i corresponding to the effective dynamic friction coef- 
ficient and tn^o corresponding to the normal stress in 
the limit of zero shear stress. This function is fit to the 
stress measurements for each viscosity curve shown in 
Figs. [2 [3l m and [51 The values of and tn^o are 
plotted in Figs. [12] [131 respectively. The error bars in- 
dicate the relative rms deviation of the linear fit from 
the data {{tn^o + /i~"^r — tn)'^)^^'^ /[max{r) — minir)]^ so 
any strongly non-linear relationship would be revealed by 
large error bars. 

The fit values of in Fig. [12] tend to be small at low 
packing fractions, then increase and level off at around 
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FIG. 13: (color online) Zero-shear rate limit of the normal 
stress tn,o from a fit of tn = rN,o + iJi~^T. Data are shown 
for several particle shapes listed in the key vs. normalized 
packing fraction. Negative stresses for < (/)c are expected 
due to surface tension at the edge of the suspension pulling the 
rheometer plates together, while positive stresses for > (/)c 
are expected due to jamming of the suspension against the 
plates. 



0.5 as (j) approaches (j)c from below, then suddenly drop 
again below 0c- The linear relationship between shear 
and normal stress, with values of the friction coefficient 
fi on the order of 1 are typical of a frictional relation- 
ship ^ommonly seen in Discontinuous Shear Thickenin- 
ing p, [23, [HI- It is notable that we never measure a 
significantly negative /i in Discontinuous Shear Thicken- 
ing, which might be expected from purely hydrodynamic 
models for shear thickening [29| . The variation of ji with 
(j)/(j)c suggests a couple of different transitions. In the 
Newtonian hmit at low packing fractions, is expected 
to go to zero as the shear flow should not produce any 
normal force in a Newtonian fluid. The increase in 
with (j)/(l)c suggests an increasing non-Newtonian contri- 
bution to the stresses up to a dense limit. A similar tran- 
sition to positive normal stresses which are linear in the 
shear stress at high packing fractions has been seen pre- 
viously for spherical particles [28^ . This transition where 
frictional forces are significant corresponds to the same 
packing fraction range of ^ ^ 0.8^c where Discontinu- 
ous Shear Thickening tends to be become dominant over 
other viscous or inertial forces as in Fig. [6] or [29]. This 
correspondence supports the idea that frictional forces 
are responsible for Discontinuous Shear Thickening [5]. 
At lower (/), it is presumably more difficult to produce 
system-spanning chains of solid contacts between parti- 
cles. 



The fit values of the static limit of the normal stress 
rAr,o in Fig- [IS] are all negative for (j) < (j)c- Observa- 
tion of the suspensions when loaded into a parallel plate 
rheometer reveals a concave curvature of the liquid-air 
interface at the side, which produces stress from surface 
tension to pull the rheometer plates together. The scale 
of this stress is expected to be about —j/d ~ —50 N, on 
the same order as seen in Fig. [131 On the other hand, in- 
creasingly positive normal stresses are seen in the static 
limit for > 0c- Above the jamming transition, particles 
are expected to push back against each other as they are 
confined to a tight packing [11]. As a result, they are also 
expected to push against the top plate, resulting in the 
positive normal forces observed. The observations of this 
transition to positive normal stresses and the change in 
friction coefficient /i at the same packing fraction where 
the yield stress jumps (our definition of 0c) are consis- 
tent with expectations that both shear and normal stress 
scalings change sharply at the jamming transition [18|. 



VI. DISCUSSION 

Here we comment on some earlier works which we may 
be able to understand better based on our more detailed 
studies of particle shape effects on the strength of Dis- 
continuous Shear Thickening. 

Clarke measured viscosity curves for several different 
shapes at a fixed packing fraction [l2[ . He reported the 
order of steepest to shallowest slope of viscosity curves: 
rods, plates, grains, spheres. Based on this he concluded 
that more anisotropic shapes shear thicken more strongly. 
In comparison, we showed in Fig. [6] that the steepness of 
viscosity curves of different convex shapes collapse when 
the packing fraction is normalized by the critical pack- 
ing fraction (j)c such that the steepness becomes greater 
approaching 0c from below. This implies that when dif- 
ferent shapes are compared at the same packing frac- 
tion, the ones with the smallest 0c will be closer to 0c 
and therefore steeper. Based on these results it seems 
more complete to identify the effect of particle shape as 
changing the value of the critical point 0c, where more 
anisotropic shapes tend to produce lower values of 0c. 

Suspensions have been found to exhibit weaker shear 
thickening with increasing particle size polydispersity at 
a fixed packing fraction [30]. In contrast, we can collapse 
data in Fig. [6] as a function of normalized packing frac- 
tion without regard to polydispersity, which ranged from 
very uniform for the PRINT® process, to very polydis- 
perse cornstarch suspensions. We suggest that the effect 
observed by Ref. [30[ was a shift in the value of 0c, which 
is known to shift upwards for more polydisperse packings 
[22]. In fact, Ref. [30] reported the value of 0c based on 
sediment concentrations, and performed measurements 
over a range of packing fractions, but still measured much 
closer to 0c for the low polydispersity suspensions than 
for the high polydispersity suspensions. Thus, the weak- 
ening observed by Ref. [30] for increasing polydispersity 
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at fixed packing fraction can be explained as a result of 
moving further away from (j)c. 

VII. CONCLUSIONS 

We measured the rheology of dense suspensions of fab- 
ricated PEG particles with different shapes, and found 
Discontinuous Shear Thickening that is qualitatively sim- 
ilar for each particle shape (Figs. [21 [3l IH and [5]). The 
logarithmic slopes of the stress/shear-rate relation in the 
shear thickening regime for convex particles are found to 
collapse onto a single curve when the packing fraction cj) 
is normalized by (pc-, which corresponds to the jamming 
transition measured independently as the onset of a large 
yield stress (Fig. [6]). Additionally, these slopes diverge in 
the limit as (j) approaches (j)c from below, as was found 
for nearly spherical particles [10], suggesting the jamming 
transition controls Discontinuous Shear Thickening like 
a critical point. 

The fact that acts as the control parameter for the 
strength of Discontinuous Shear Thickening has several 
consequences. It suggests that the apparent weakening 
of shear thickening for more spherical shapes [12] or more 
polydisperse sizes [30| at constant (j) can be explained by 
the fact that (j)c is moved further away. Similarly, par- 
ticles that form fractal aggregates have been found to 
exhibit Discontinuous Shear Thickening at packing frac- 
tions as low as 11%, which is not surprising in this in- 
terpretation since they also pack at very low densities 
(31j . For designer suspensions, this implies that strong 
shear thickening can be achieved using less material by 
designing particles that pack at low densities. 

For one shape we studied, the hooked rods, the slopes 
of the viscosity curves did not collapse with the same 
scaling as for convex particles. We can only speculate 
on why hooked rods show a different scaling. The non- 
convex shape allows for more contact points, which can 
frustrate motions in both shear and compressional direc- 
tions on contact, rather than just compression for convex 



particles. The hooks may frustrate particle rotation and 
passing in tight packings under shear flow, but it is not 
clear how these features specifically translate to the wider 
packing fraction range for strong shear thickening. 

We identified two other mechanical properties which 
relate Discontinuous Shear Thickening to jammed sys- 
tems. First, the shear stress and normal stress are lin- 
early related as in a frictional scaling for (j) > 0.8(/)c 
where Discontinuous Shear Thickening and jamming oc- 
cur (Figs.[TTl[T2]), which suggests that energy dissipation 
comes from friction between solid particle contacts that 
span the system in both cases, rather than from viscous 
dissipation. Additionally, we showed the scale of the yield 
stress comes from interfacial tension for jammed suspen- 
sions, which can be seen directly by observing the surface 
roughness of the suspension (Fig. [8l [9l [TQ|) . This is sim- 
ilar to when the suspensions dilate against the liquid-air 
interface to result in Discontinuous Shear Thickening, al- 
though generally the forces need not come from interfa- 
cial tension, but could come from any forces which con- 
fine the system at the boundary |5|]. The difference be- 
tween the jammed state above (j)c and the Discontinuous 
Shear Thickening regime is that the latter requires dila- 
tion from shear to achieve the solid particle contacts to 
reach the frictional state, while the jammed state already 
has solid particle contacts at rest. The similarity in the 
mechanics of force transmission suggests the possibility 
that a generalization of jamming theory may also be able 
to describe Discontinuous Shear Thickening systems, as 
was suggested by [9.]. 
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